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Abstract 

The /3-decay process of the n Li halo nucleus into 9 Li and d is studied in a three-body model. The 
n Li nucleus is described as a 9 Li + n + n system in hyper spherical coordinates on a Lagrange mesh. 
Various 9 Li + d potentials are compared involving a forbidden state, a physical bound state, and 
a resonance near 0.25 MeV in the s wave. With an added surface absorption, they are compatible 
with elastic scattering data. The transition probability per time unit is quite sensitive to the 
location of the resonance. For a fixed resonance location, it does not depend much on the potential 
choice at variance with the 6 He delayed deuteron decay. The calculated transition probability per 
time unit is larger than the experimental value but the difference can be explained by a slightly 
higher resonance location and/or by absorption from the 9 Li + d final channel. 
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I. INTRODUCTION 



The (3 decay with emission of a deuteron, also known as (3 delayed deuteron decay, was 
first observed for the 6 He halo nucleus [H, Q, [3|. The difficulty of the measurement led to 
conflicting results and raised a number of theoretical questions. The problems are related 
with the fact that the branching ratio is much smaller than expected from simple i?-matrix 
[l| or two-body 0] model estimates. 

A semi-microscopic study of the process has been able to explain that the low value of the 
branching ratio is the result of a cancellation between contributions of the "internal" and 
"halo" parts of the Gamow- Teller matrix element A fully microscopic description of the 
f3 decay of the 6 He nucleus to the a + d continuum 0] then indicated that the cancellation 
could be more important than expected from the available data of that time. Both works 
emphasized that the (3 delayed deuteron decay probes the halo up to rather large distances. 

These interpretations were confirmed in a recent study based on a three-body model 0] ■ 
It was shown that the cancellation requires that the a + d potential contains a forbidden 
state below the 6 Li ground state in order to have the correct node structure of the scattering 
wave function. The cancellation is so strong in the Gamow- Teller matrix element for 6 He 
that it requires an almost perfect balance between the internal and halo parts, which should 
be fortuitous. A similar effect is thus not expected for other halo nuclei possessing a (3 
delayed deuteron decay branch such as n Li or 14 Be. 

The most interesting halo nucleus is probably n Li j^]. Its two-neutron separation energy 
is particularly small: 300 ± 19 keV according to the atomic mass evaluation [9j or 376 ± 5 
keV according to a recent preliminary result [10]. This nucleus can be considered as a 9 Li 



core surrounded by two halo nucleons distant from each other by more than 6 fm [ll|j m 
agreement with theoretical expectations [l2| . It differs from 6 He by the fact that its core 
does not correspond to a closed shell. The halo structure is understood as due to a virtual 
state in the s wave of the 9 Li+n interaction [l3|]. Another difference with 6 He is that the core 
is unstable. Therefore many more (3 decay channels are o pen . This complicates experiments 
but also offers many opportunities to test models [H, 15. lid]. Among the possible channels, 
the delayed deuteron (3 decay remains especially interesting because this decay essentially 
occurs inside the halo and can probe its properties. In experiments however, the deuteron 



decays can not easily be separated from the delayed triton channel [17j]. Anyway the (3 
delayed deuteron decay has been observed lU, 16[ with a branching ratio of (1.5 ±0.2) x 10~ 4 
16[ . This order of magnitude is consistent with predictions of a simple model [l8[ and of a 



model based on a limited hyperspherical-harmonics expansion [19j. The analysis of a new 



experiment is in progress [20J . It is thus timely to reexamine this process at the light of the 
recent knowledge gained on the 6 He decay. 

In the present work, we describe 11 Li as a 9 Li + n + n system with effective 9 Li + n and 



n+n interactions treated in hyperspherical coordinates with the Lagrange-mesh method [21 
The final 9 Li + d state is modeled as a deuteron wave function multiplied by a scattering 
wave function deduced from a potential. Experiments have revealed the important role 
played by a resonance around 18 MeV in the excitation spectrum of n Be [l6|. Taking also 
that information into account, the transition probability per time and energy units will be 
calculated. The total transition probability is constrained with the branching ratio. They 
will be analyzed by comparison with the 6 He decay with emphasis on the role of the node 
structure of the scattering wave functions. 

In Sec. [Tll the /3-decay model for the n Li two-neutron halo nucleus into the 9 Li + d 
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continuum is summarized. The potentials and the corresponding three-body hyperspherical 
and two-body scattering wave functions are described in Sec. IIHI The properties of the 
Gamow- Teller matrix element are studied in Sec. [TV] In Sec. [V] we discuss the obtained 
numerical results and compare them with experimental information. Conclusions are given 
in Sec. EQ 



II. MODEL 

The initial wave function is expressed as a bound-state wave function of the three-body 
9 Li + n + n system with local 9 Li + n and nn interactions. We neglect the spin of the core. 
The total orbital momentum L of the three particles is assumed to be equal to the total spin 
S of the neutrons as for 6 He 0. Jacobi coordinates, i.e. the relative coordinate r between 
the neutrons and the coordinate R of their center of mass with respect to the 9 Li core, 
are necessary to calculate the overlap with the final scattering state. These coordinates 
are conveniently replaced by hyperspherical coordinates which involve five angular variables 
^5 and the hyperradius p. The wave function is expanded over hyperspherical harmonics 
depending on Q 5 and on the hypermomentum K. The coefficients in this expansion depend 
on the hyperradial coordinate p and are expanded in Lagrange functions j22| (see Ref. [2l[ 
for details). 

For the final scattering state, we assume an expression factorized into the deuteron 
ground-state wave function depending on r and a 9 Li + <i scattering wave function depending 
on R derived from a potential model. We neglect the small D component of the deuteron. 

The transition probability per time and energy units is given by |23j 

dW m P c 2 



lis « 

where m e is the electron mass, v and E are the relative velocity and energy in the center 
of mass system of 9 Li and deuteron, and Gp = 2.996 x 10 -12 is the dimensionless /?-decay 



constant [24(. The Fermi integral f(Q — E) depends on the kinetic energy Q — E, available 
for the electron and antineutrino. The mass difference Q between initial and final particles 
is given in MeV by 

Q = 3.007 - S 2n (2) 

as a function of the two-neutron separation energy of the halo nucleus. With the n Li 
value S2n = 300 ± 19 keV from the atomic mass evaluation [9], Q is equal to 2.71 ± 0.02 
MeV. However according to a recent remeasurement, the n Li two-neutron separation energy 



becomes S 2n = 376 ± 5 keV [lOj leading to Q = 2.63 MeV. We shall first use the standard 
value and then consider the importance of this modification. 

Since the total orbital momentum and parity are conserved, only the I = partial scat- 
tering wave contributes. Hence, only the initial L = S = component of n Li described with 
a spin core can decay to 9 Li + d. In order to allow the use of a complex optical potential 
for describing the scattering states, we generalize the formula of Refs. p, 0] • The final state 
is described by an ingoing scattering wave. At energy E, a partial wave u~ E \ of an ingoing 

scattering wave function is related to a partial wave u^J of an outgoing scattering wave 
function by 

u { e}{R) = (-1)'41*(#)- (3) 
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The outgoing radial scattering wave functions 

u { +}(R) = e iSl u E>l (R) (4) 
are normalized asymptotically according to 

u E ,i{R) -> cos^(^)^(^) + sin^(^)^(^), (5) 

R— >oo 



where k is the wave number of the relative motion, Fi and Gi are Coulomb functions [25 
and 5i(E) is the /-wave phase shift at energy E. The subscript / = is understood in the 
following. 

The reduced transition probability can be written as 



B GT (E)=6\ 2 \e l5 °I E (oo)\ 2 , (6) 



where A = —1.25 26[. The phase in front of I E does not play any role if the potential is 
real. The integral 

I E (R) = / u E {R')u eS {R')dR' (7) 
Jo 

depends on a cutoff radius R over the relative coordinate between the core and the center 
of mass of the nucleons. Only its value at infinity is physically relevant but it will help us to 
understand the physics of the decay process. This integral involves scattering wave functions 
u E (R) and depends thus on the 9 Li + d relative energy E. This integral also involves an 
effective wave function 

poo 

u cS (R)=Rj2 Z K (r,R)u d (r)rdr, (8) 
k Jo 

where Ud(r) is the deuteron radial wave function depending on the relative coordinate r of 
the two nucleons. The sum runs over the values of the hypermomentum K in the expansion 
of the initial bound state. The function Zk(t,R) is the radial part of the K component 
with all angular momenta equal to zero in the expansion of the initial wave function. Its 
expression is given by Eqs. (3) and (14) in Ref. M]) where however a normalization factor 
[(A — 2)/A] 3//4 is missing. The results of Ref. [7J) must be modified accordingly. In the 
following, we also make use of partial integrals /^(i?) obtained from Eq. ([7]) with the 
different terms in Eq. (jSJ). The sum of the lf\R) is I E (R). 



III. POTENTIALS 



The deuteron wave function Ud is calculated with the central Minnesota interaction 27 
'see Ref. |2l| for details). An energy Ed = —2.202 MeV is obtained. 



The 9 Li+n+n wave function is calculated with the 9 Li+n potential P2 of Ref. 13j and the 
nn Minnesota interaction with exchange parameter u = 1. In order to fit the binding energy 
of 11 Li, the P2 interaction is multiplied by a parameter [121 ] . The values 0.992 and 0.9965 
provide S2 n = 0.307 and 0.376 MeV, respectively. The s-wave scattering length is then 
slightly modified from —25.4 fm to —19.0 or —22.2 fm, respectively. Potential P2 contains 



4 



a forbidden state in the s wave which is eliminated with the pseudo-potential method 28 . 
Forbidden states need not be eliminated in two-body systems as they do not affect scattering 
properties. Their presence leads to more realistic wave functions for the relative motion. In 
three-body systems however, forbidden states must be eliminated because otherwise they 
would unrealistically contribute to the binding energy. The pseudopotential moves them to 
a high energy without affecting the other properties of the two-body potentials. 

A 9 Li + d optical potential has been obtained by fitting elastic scattering data at a cm. 
energy of 3.86 MeV j29|. The real part of this potential does not display any resonance 
below the Coulomb barrier. Such a resonance has been observed in several channels at the 
excitation energy 18.15 ± 0.15 MeV i.e. at the cm. energy 0.25 ± 0.15 MeV above 



the 9 Li + d threshold. As shown below, this resonance is crucial to explain the order of 
magnitude of the (3 delayed deuteron decay of 11 Li. The potential of Ref. [29] is thus not 
useful here. Its real part provides three bound states. When the depth of its real part is 
reduced from 104.6 MeV to 89 MeV, the upper bound state becomes a resonance near the 
experimental value. However the agreement with the elastic scattering experiment is then 
lost. 

We approximate the 9 Li + d potential by expressions based on simple physical arguments 
derived from a microscopic cluster model interpretation, (i) At short distances, 9 Li and 
deuteron can form a bound state in the s wave. This bound state has the same parity as 
the 9 Li core, i.e. a negative parity. We thus impose to the potential to reproduce the energy 
of the l/2~ excited state of n Be at an excitation energy of 0.320 MeV. This means that our 
potentials will have a bound state near —17.6 MeV. (ii) In the microscopic cluster model, 
the 9 Li + d system possesses a forbidden state in the s wave. The role of such a state can 
be simulated by a potential deep enough to contain an unphysical bound state below the 
physical bound state in order to simulate the correct node structure of the scattering wave 
function, (iii) The n Be nucleus displays a resonance around 0.25 MeV above the 9 Li + d 
threshold. 

As in Ref. jl8], we consider simple Gaussian potentials parametrized as 

V(R) = -V exp(-aR 2 ). (9) 

The choice of a Gaussian form factors restricts the number of parameters. A Coulomb term 
3e 2 erf (/3r)/r with (3 = 0.71 fm _1 (scaled from 0.75 fm _1 in the a + d case) is added to all 
potentials. The Coulomb barrier is located between 0.55 and 0.6 MeV. 

Potential V a with a = 0.14 fm -2 and Vq = 89.5 MeV (see Fig. [T]) has a bound state 
at energy —17.63 MeV and a forbidden state at —52.42 MeV. The width of the Gaussian 
form factor has been chosen in such a way that the potential also verifies criterion (iii). The 
corresponding phase shift is displayed as a full line in Fig. [2j A resonance appears at about 
0.33 MeV with a width of about 0.1 MeV. The spin and parity of this resonance should be 
3/2~. Its width is smaller than the experimental width derived in Ref. [16j which however 
largely exceeds the Wigner limit and is therefore questionable. 

In the 6 He forbidden state plays a crucial role in the reproduction of the ex- 

perimental order of magnitude. In order to study the role of the forbidden state here, we 
perform pairs of supersymmetric transformations [30[ in order to remove it from V a while 
keeping the other bound state and the s-wave resonance and phase shift. The resulting 
phase-equivalent potential denoted as V a i exhibits a strong repulsive core (see Fig. [T]). The 
physical bound state of V a \ is then removed by another pair of transformations leading to 
the phase-equivalent potential V a 2 without any bound state. Both potentials V a \ and V a 2 
provide the same s-wave phase shift as V a in Fig. [2j 
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FIG. 1: Potential V a and phase-equivalent potentials V a \ and V a 2 (full lines); potential Vj, (dashed 
line); potential V c and phase-equivalent potentials V c i (dotted lines). 

We also consider other Gaussian potentials. Potential V& with the same range as V a but 
Vo = 42.7 MeV has its ground state at —17.63 MeV. This potential possesses a weakly 
bound state near —0.184 MeV in place of a resonance. Potential V c is quite similar to V& 
but differs from it by the fact that it possesses a resonance at 0.28 MeV in addition to a 
bound state at —17.66 MeV, with a = 0.161 MeV and V = 44.8 MeV. Removing the bound 
state leads to the phase-equivalent potential V& with a repulsive core. These potentials are 
compared in Fig. [T] and their phase shifts are displayed in Fig. [2j One observes that the 
phase shift of potential V c has the same shape as the phase shift of V a and also displays the 
expected resonance but at a lower energy. On the contrary, the phase shift of potential V& 
is monotonic. 



These potentials are compatible with the elastic scattering data of Ref. 2j| if some surface 
absorption is added. Without absorption, even the order of magnitude of the cross section 
is incorrect beyond 70 degrees. We use the simple optical potential 

V opt (R) = -(V + iW ^R)e- aR2 , (10) 

where the imaginary part is proportional to the derivative of the real part. The ratio of the 
elastic cross section to the Rutherford cross section is compared with experiment in Fig. [3] 
for an optical potential based on the parameters of V a . The only free parameter is thus Wq. 
One observes that a fair agreement with the data is obtained for W = 50 ± 10 MeV, except 
around the first maximum near 40 degrees. The present real part is thus compatible with 
elastic scattering. 
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FIG. 3: Ratio of the elastic cross section to the Rutherford cross section at a cm. energy of 3.86 
MeV for potential (HQ) with a = 0.14 fm~ 2 , V = 89.5 MeV, and Wn = 40 (dashed line), 50 (full 
line), and 60 (dotted line) MeV. Experimental data are from Ref. [291 ]. 

IV. GAMOW-TELLER INTEGRALS 

The integrals 1^ (R) calculated with potential V a are displayed in Fig. @] as a function 

of R for different K values at energy 1 MeV. The convergence of YIk ^e^ 1 (-^) ^ s reached for 
-Kmax = 20. Partial waves K = 2, 0, and 4 are strongly dominant, although the cumulated 
contribution of all higher partial waves is not negligible. Contrary to the 6 He case, no 
important cancellation is encountered here because the dominant contributions have the 
same sign. 

In Fig. [51 the integrals Ie{R) are compared for the different potentials considered at two 
energies: the near-resonance energy 0.3 MeV and a typical non-resonant energy 1 MeV. As 
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FIG. 4: Partial integrals 1^ i R ) for K = 0, 2, and 4 at the energy E = 1 MeV for potential 
V a . The sum of the three components K = 0, 2, and 4, (dotted line) and the converged sum 
I E (R) = Y,k i{ e\ r ) [ E( l- ©] (lowest full line) are also displayed. 



shown by Eq. (J7|), the integral displays a minimum or a maximum every time either ue 
or itcff vanishes. Its behavior depends on the node structure of the scattering wave function 
and thus on the depth of the potential. At both energies, all curves present an extremum 
near 2 fm which corresponds to the unique node of u e g. They present one to three additional 
nodes corresponding to the possible bound states and resonance of the potential. However, 
in spite of their different numbers of bound states and thus of nodes, potentials V a , V a \, and 
V a 2 do not give very different results at both energies. In all three cases, the amplitude of 
the integral starts to increase beyond 3 fm and reaches a plateau near 20 fm. At 1 MeV, 
it presents a maximum near 9 fm. This maximum is at the same location for V a , V a \, and 
V a 2 because phase-equivalent potentials have the same asymptotic behavior and thus the 
same nodes beyond the potential range. This situation must be contrasted with the 6 He 
case where the cancellation enhances tiny differences and where phase-equivalent potentials 
provide very different results 0]. 

The results for potential V c which satisfies the same physical conditions as V a \ are very 
similar because the scattering wave function has the same number of nodes and similar 
locations of these nodes. At 0.3 MeV, the integrals have opposite signs for V a and V c beyond 
4 fm because the resonance is below 0.3 MeV for V c while it is above for V a . They have the 
same sign and similar magnitudes at 1 MeV. 

On the contrary, the results obtained with V& are very different, even off resonance, 
because the scattering wave function has nodes at quite different locations. In particular its 
node near 8 fm at 0.3 MeV or 6 fm at 1 MeV leads to a cancellation similar to that of the 
6 He case. We shall see in the next section that this type of result is ruled out by experiment. 
It is important to realize that has the same physical bound state as V a \ and V c near —17.6 
MeV. However, V h does not reproduce the resonance (see Fig. [2]). 



V. TRANSITION PROBABILITY PER TIME AND ENERGY UNITS 



The transition probability per time and energy units given by Eq. ([T]) is plotted in Fig. E] 
as a function of the relative 9 Li + d energy E for the different potentials. For the potentials 
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FIG. 5: Integrals I E (R) [Eq. ©] at (a) E = 0.3 MeV and (b) 1 MeV, for various potentials. 



displaying a resonance, the results are qualitatively very similar. The shape of the curve 
does not depend much on the resonance location. On the contrary, potential V& provides 
results with a much smaller order of magnitude. The same situation is observed with other 
potentials that do not possess a resonance. 

The total transition probabilities per time unit (integrated from or from some cutoff 
to Q) corresponding to the various potentials are compared with the experimental value in 
Table [H This value is calculated from the experimental branching ratio 1Z by 

W exp = ^ln2/t 1/2 Ri81.57es- 1 , (11) 

where £x/2 is the n Li ground-state half life 8.5 ms. The experimental branching ratio is 
(1.5 ± 0.2) x 10~ 4 [16]. The results at various cutoff values in Table fl] will be useful for 



comparison with the new experiment [20 



All potentials except V& provide the right order of magnitude but overestimate the exper- 



imental value of Ref. [16j by a factor larger than 3. Results of successive phase-equivalent 
potentials differ by about 10 %. The main difference between V a and V c arises in the cutoff 
dependence which is sensitive to the resonance location. 

Table [J also indicates the role of a larger two-neutron separation energy S^. This in- 
troduces a modification of the Q value and of the n Li wave function. Except for the 
transition probabilities are slightly reduced, by about 20 %. This effect is rather weak and 
does not modify the discussion. The V& variation emphasizes the high sensitivity to weak 
modifications when a cancellation occurs, like in the 6 He case. 
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E (MeV) 



FIG. 6: Transition probability per time and energy units dW/dE of the n Li decay into the 
9 Li + d continuum as a function of the relative 9 Li + d energy E calculated with various 9 Li + d 
potentials. 

TABLE I: Total transition probability per second W (in 10~ 3 s _1 ) for the (5 decay of n Li into 
9 Li + d. For each value of the two-neutron separation energy Si n (in MeV), the rows correspond 
to various cutoffs. 

S 2n cutoff V a Vgi V a2 V b V c Vd Exp. 
0.307 E > 38.1 42.1 46.7 0.0718 59.7 54.1 12 ± 2 [16] 

E > 0.3 31.0 34.3 38.2 0.0392 22.0 19.7 

E > 0.5 4.7 5.4 6.2 0.0096 2.7 2.3 
0.376 E > 31.5 34.9 39.0 0.1014 50.8 45.8 12 ± 2 [16] 

E > 0.3 25.7 28.6 32.1 0.0622 19.0 17.0 

E > 0.5 4.0 4.6 5.3 0.0185 2.4 2.1 



Let us study the role of the main uncertainties in our theoretical description. The first 
uncertainty concerns the energy location of the resonance. The location of the peak in 
Fig. affects the total transition probability. In Table [TTJ we study the dependence of the 
transition probability on the resonance energy E r . To this end, we slightly vary the depth 
Vo in potential V a . This leads to a small violation of our criterion (i), i.e. the energy Ebs 
of the physical bound state is somewhat modified, but this modification remains acceptable 
in view of our other simplifying assumptions. As shown by Table [Til W is locally quite 
sensitive to the resonance energy and a slightly higher location would lead to smaller values. 
A higher location of the resonance also reduces the cutoff dependence. 

Another effect, not encountered in the 6 He case, arises from the fact that several channels 
are open below the 9 Li + d channel, the lowest one being the 10 Be + n channel. Transfer 
towards these channels is possible at all energies but should be rather weak below the 
Coulomb barrier. The magnitude of the surface absorption into these channels has been 
derived from experimental scattering data in Fig. [3] but at an energy 3.86 MeV much higher 
than the Coulomb barrier. Therefore we restrict the discussion to small values of Wq in 
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TABLE II: Dependence of the total transition probability per second W (in 10~ 3 s _1 ) on the 
resonance energy E r (in MeV) calculated with a Gaussian potential with a = 0.14 fm~ 2 as a 
function of its depth Vo (in MeV) for various cutoffs. The forbidden state energy Eps an d the 
physical bound state energy Ebs are also displayed. 







Ebs E r 




W 










E > E > 0.3 E > 0.5 


90.8 


-53.44 


-18.26 0.25 


47.5 


10.7 


2.1 


90.1 


-52.89 


-17.92 0.30 


42.1 


24.6 


3.3 


89.5 


-52.42 


-17.63 0.35 


38.1 


31.0 


4.7 


89.0 


-52.02 


-17.38 0.40 


35.2 


31.0 


6.2 


88.5 


-51.63 


-17.15 0.46 


32.5 


29.6 


8.0 



Eq. ffTUl) (this parameter should probably depend on energy but we neglect this effect here). 
One observes in Fig. [7] that the role of the resonance is strongly reduced even by a weak 
absorption. On the contrary, the results above 1 MeV are not much affected. The energy 
dependence of dW/dE becomes weaker when absorption increases. 

i 



a; 1(T 2 

1 

^ io- 4 
i(r 6 

0.5 1 1.5 2 2.5 

E (MeV) 

FIG. 7: Transition probability per time and energy units dW/dE of the 11 Li f3 decay into the 
9 Li + d continuum as a function of the relative 9 Li + d energy E calculated with various values of 
the surface absorption strength Wq (in MeV) in Eq. (flQ|) . 

The total transition probabilities per second W calculated with potential ( fTUl) in Eq. ([6]) 
are displayed in Table III II with a and Vq as in potential V a for several fixed values of the 
surface absorption constant Wq. One observes that a much smaller absorption than in Fig. [3] 
leads to a strong reduction of the transition probability. As explained by Fig. [71 absorption 
leads to a weaker dependence on the cutoff. 

VI. CONCLUSIONS 

In the present work, we studied the /9-decay process of the n Li halo nucleus into the 9 Li+d 
continuum in the framework of a three-body model. Three-body hyperspherical bound-state 
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TABLE III: Dependence of the total transition probability per second W (in 1CP 3 s _1 ) on the 
surface absorption strength Wo (in MeV), calculated with potential V a as real part for various 
cutoffs. 

Wo E > E > 0.3 E > 0.5 






38.1 


31.0 


4.7 


1 


16.7 


13.0 


3.5 


2 


9.9 


7.8 


2.7 


5 


4.4 


3.6 


1.8 



wave functions on a Lagrange mesh and two-body 9 Li + d scattering wave functions have 
been used. For the calculation of the /3-decay transition probabilities per time and energy 
units, several 9 Li + d potentials were employed. 

Some 9 Li+<i potentials are physically inspired by a microscopic cluster picture and involve 
a forbidden state and a physical bound state simulating the l/2~ excited state of 11 Be. A 
resonance occurs in the s wave at about the experimental energy. For potentials of this 
family, the transition probability per time unit is weakly sensitive to the potential choice. 
However a potential without this resonance fails to reproduce even the order of magnitude 
of the transition probability. The high sensitivity of the delayed ft decay of 6 He due to a 
cancellation in the Gamow- Teller matrix element does not occur here. This is emphasized 
by using phase-equivalent potentials differing by their number of bound states: they give 
very different results for 6 He and very similar results for n Li. 

The theoretical result is sensitive to the 11 Li separation energy (about 20 % if S 2n is 
increased by about 70 keV). It is also sensitive to the location of the resonance. A more 
accurate experimental determination of the location of this resonance (with a precise def- 
inition of the resonance energy for such a broad resonance) would be very useful. Elastic 
scattering data on d( 9 Li, 9 Li)<i at an energy close to the Coulomb barrier extending up to 
backward angles might also help reducing the uncertainty on the parameters of the optical 
potential for this collision. 

The calculated transition probability overestimates the experimental result of Ref. [l6[ 
by a factor larger than 3. The overestimation can be reduced by modifying the resonance 
location and/or by introducing absorption removing flux from the 9 Li + d final channel. 
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